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Abstract. We present the latest version 2.6 of FeynHiggs, a program for computing MSSM Higgs- 
boson masses and related observables, such as mixing angles, branching ratios, and couplings, 
including state-of-the-art higher-order contributions. The most important new feature is the inclu- 
sion of the fully complex 0{atas) two-loop corrections, which enables FeynHiggs to give the most 
precise Higgs-mass evaluation in the complex MSSM in the Feynman-diagrammatic approach to 
date. 

PACS. 12.60.Jv Supersymmetric models - 14.80.Cp Non-standard-model Higgs bosons 



1 Complex Parameters in the MSSM 
Higgs Sector 

The Higgs sector of the Minimal Supersymmetric Stan- 
dard Model with complex parameters (cMSSM) con- 
sists of two Higgs doublets 



H, = 



^(<^i-iXi) 



(1) 



(2) 



+ 71 ((^2 +1X2) 

which form the following Higgs potential 



V = ml HiHi + ml H^H-z - {m\^ e^i^H'^B^ + h.c.) 

^^l±A^U,Hi-R2H2f + ^^\R,H2\\ (3) 
o z 

The Higgs potential contains two complex phases ^, 
arg(mf2)- The phase &Y%{vfi\^ can be rotated away 
[1115] and, at tree level, has to vanish In order to 
fulfill the minimum condition of the Higgs potential, so 
there Is no CP-vlolatlon at tree level and the spectrum 
contains five states of definite CP-parity: /i, iJ, A, H^. 
In the following we review the Inclusion of higher-order 
corrections to Higgs-boson masses and more Into the 
code FeynHiggs J3i|4l[5l[6] . 

CP-vlolatIng effects are Induced by complex pa- 
rameters that enter via loop corrections: the Hlggslno 
mass parameter fi, the trlllncar couplings At,b^T-, and 
the gaugino mass parameters Mi_2.3. They yield ShA, 
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Sha 7^ and Induce mixing between h, H, and A [7]. 
The Higgs mass matrix has the form 



2 

2^Hh q 



Sha 



(4) 



Sai 



— m\ + ^AA J 



where mh^H.A denote the tree-level Higgs masses, and 
it should be noted that in general is symmetric 
but not Hermltian. In the approximation of vanish- 
ing external momentum (q^ = 0), one can obtain the 
higher-order corrected mass eigenstates via a unitary 
transformation from the tree-level states: 



'hi\ Un Uv2 Ui3 

/l2 = U21 U22 U23 

. hsi \U3i U32 U33 




(5) 



2 Higgs-boson self-energy corrections in 
FeynHiggs 

2.1 Higgs-boson masses 

The following contributions to the mass matrix and 
the charged-Hlggs-boson self-energy are taken into ac- 
count: 
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® Leading 0{atas) cMSSM two-loop corrections [8]. 
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® Leading 0{af) and subleading 0{abas, atab, a^) 
two-loop corrections evaluated in the MSSM with 
real parameters (rMSSM), where the phases are 
included only partially ^[TOllllj. 

® Full one-loop evaluation (all phases, dependence) 
[6j and leading non-minimal flavour-violating 
(NMFV) corrections [12] . 

FeynHiggs performs a numerical search for the com- 
plex roots of detM^((7^) which are denoted as M^., 
i — I . . .3. A decomposition can be performed, 

M'^ ^ -iMT, (7) 

where M is the mass of the particle and F its width. 
We then define the loop-corrected masses according to 

Mh, ^ Mh, s=; Mh, . (8) 

The Higgs masses are thus determined as the real parts 
of the complex poles of the propagator. Complex con- 
tributions to the Higgs mass matrix (from Im S) are 
taken into account ^6, 13 . The diagonalization routines 
are available as a stand-alone package from the Web 
site www.feynarts.de/diag [M]. 

2.2 Two-loop corrections in the complex MSSM 

Including the phase dependence, the complete one- 
loop [6] and the two- loop contribution of 0{atas) [8j to 
the Higgs self-energies are taken into account. Within 
the Higgs sector, the parameters have to be defined up 
to 0{atas)- The masses of the charged Higgs boson, 
the Z-boson, as well as the W-hoson are defined as 
pole masses, 

6M\^'^ = ReZ'jj|Y(M|) with X = {H^ , W, Z}, (9) 

with [i) denoting the loop order. Furthermore, it is 
required that there be no shift of the minimum of the 
Higgs potential which is fixing the tadpole parameters, 

5tf = -rf with (f) = {/i, H, A] . (10) 

The Z- factors and tan/3 are defined within the DR- 
scheme [TS l lTB ] . 

The parameters of the top sector have to be defined 
at one-loop level. The top-quark mass and the top- 
squark masses are fixed by an on-shell condition and 
the mixing angle and the corresponding phase by 

Re (m? J -f Ite (m| ) = , (11) 

generalizing the renormalization conditions imposed in 
|17j for the use of complex parameters. 

To extract the relevant terms at two-loop order we 
used the approximation of vanishing external momenta 
and vanishing electroweak gauge couplings in the eval- 
uation of all two-loop diagrams including those needed 
for calculating the two-loop counterterms. 

A new flag in FeynHiggs (see Sec. 16.21 below) con- 
trols the treatment of phases in the part of the two- 
loop corrections known only in the rMSSM so far. The 
following options are possible: 



— all corrections: ©(ata^, ctbas, , atCKfj, a^) in the 
rMSSM, 

— only the cMSSM 0{atas) corrections, 

— the cMSSM 0{atas) corrections combined with 
the remaining corrections in the rMSSM, truncated 
in the phases, 

— the cMSSM 0{atas) corrections combined with 
the remaining corrections in the rMSSM, interpo- 
lated in the phases [default]. 

FeynHiggs thus not only has the most precise eval- 
uation of the Higgs masses in the cMSSM available 
to date (using the Feynman-diagrammatic approach), 
but also a method to obtain a reasonably objective 
estimate of the uncertainties due to the rMSSM-only 
parts. 

Implementing the 0{atas) cMSSM corrections in 
FeynHiggs was a major piece of work. The amplitudes 
could be shrunk from 38 MB to less than 1.5 MB, 
mainly by abbreviationing techniques and exploiting 
the unitarity of the sfermion mixing matrices. The 
compile time is about 3 min (up from 45 sec in Feyn- 
Higgs 2.5) and the run time is 28 msec per parameter 
point (up from 27 msec in FeynHiggs 2.5). These fig- 
ures show that the full cMSSM evaluation is actually 
usable in everyday life. 

As a numerical example we show in Fig. [l]the de- 
pendence of Af/ij on the phase of the top-squark mix- 
ing, Xt — At — cot p. The plot shows the one-loop- 
corrected Higgs mass Mh^ as dotted curve. The drawn- 
through curve depicts the Higgs mass Mh^ including 
contributions of 0{atas)- The boundaries of the bands 
are calculated in the following way: 

M,^-(^xJ = MrrifxJ + AMnAvx. = 0) , (12) 
M^f (^x. ) = M,™"- (^x. ) + AAh, (^x. = TT) , (13) 
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0{a) : tan/? = 5 • • tan,(3 = 15 

0(a + ciiQ',) : tan /3 = 5 tan /3 = 15 

0{a + Q';a^) + interpol. : tan/^ = 5 ' ' tan/3 = 15 

{lower curves) (upper curves) 

Fig. 1. Tlie liglitest cMSSM Higgs-boson mass as a func- 
tion of tpxt at the one- and two- loop level (see text). The 
other parameters are: Msusy = M3. — M2 — 500 GeV, 
Ml = 250 GeV, fj. = 1000 GeV, M^± = 150 GeV, 
l^tl = 700 GeV. 
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where Mj^^ and M^^ respectively give the lower and 
the upper boundary of the bands for 



AMh, (ifix, - 0) ^ AMh, {(fix, - tt) . 



(14) 



If, unlike in our numerical example, Eq. (jl4p does 
not hold, Mjj°^ and AI^^ have to be interchanged in 
Eqs. and in]). 

M^°"- are the values for Af^j^ including the full one- 
loop and the 0{atas) corrections with the full phase 
dependence. AMh^ gives the size of the contributions 
that are only known for real parameters, namely those 
of 0{af,abas,atab, al). 

The crossed curve shows Mh^ , taking into account 
the 0{atas) contributions and interpolating AM^^, 
i.e. the corrections of 0{a^,abas,atCtbTOL^)- The fact 
that these crossed curves lie between the lower and 
the upper boundaries of the corresponding band shows 
that the interpolation procedure is working well. 

For the parameters chosen here (especially due to 
a relatively small value of Mfj±) the 0{atas) con- 
tributions decrease the phase dependence. As known 
from the rMSSM [i], they cause a shift of M/^ to- 
wards lower values with respect to the one-loop cor- 
rected mass, the ©(a^ , a^as, ata^, a^) corrections in- 
crease again the size of the M^^ . 



2.3 Mixing of the Higgs bosons 

FeynHiggs returns two different 'mixing' matrices. 

— UHiggs is a 'true' mixing matrix in the sense of 
being unitary and hence preserving probabilities. 
When applying effective couplings for internal Higgs 
bosons, this matrix must be used. 

It should be noted that to obtain a unitary ma- 
trix, it is mathematically a necessity that has 
no imaginary parts ~ making it Hermitian. This of 
course constrains the achievable quality of approx- 
imation. 

— ZHiggs is a matrix of Z-factors. It guarantees on- 
shell properties for external Higgs bosons [B], see 
Eq. dlHl) below. 



It is important to understand that ZHiggs and UHiggs 
are two objects with physically and mathematically 
distinct properties. Neither is universally 'better' than 
the other. 

UHiggs can be computed in two approximations: 



on-shell: Si 



;))• 



— = (see Eq. ([5|)). UHiggs coincides with ZHiggs 
in this limit and corresponds to the effective po- 
tential approach. In the absence of CP-violating 
effects, i.e. 2x2 mixing only, this is identical to 
the aeff description [T5] . 

ZHiggs is engineered to deliver the correct on-shell 
properties of an external Higgs boson, but is not nec- 
essarily unitary The following picture shows the 
type of mixing contributions which appear in the de- 
cay of an external Higgs boson (the contributions from 




mixing with the Goldstone boson and with the longitu- 
dinal component of the Z-boson are numerically small 
and hence neglected): 



h,H,A 



Disregarding possible permutations for the moment 
(but see below) , the corresponding combination of am- 
plitudes are 

= \fZfi[rh + Zhu^H + ZhA^A) (15) 

= \fZ^{ZHhrh + Th + ZnArA) (16) 
A3 = ^{ZAhFh + ZahFh + Fa) (17) 

where 

— rh,H,A is the amplitude for h,H,A-^ X, 

^ \/ ZhM,A sets the residuum of the external Higgs 
bosons to 1, 

— ZhH, ZfiA describe the transition h ^ A, etc. 

For convenience, the Z factors can be arranged in ma- 
trix form: 



ZHiggs = 



In this guise, ZHiggs can be used very much like UHiggs 
even though its theoretical origin is quite different. Re- 
assuringly, ZHiggs and UHiggs coincide in the limit 
g2 = 0. 

The transition factors Zij involve both the tree- 
level mass nii and the loop-corrected mass Aii of each 
Higgs boson: 




Zi- 



= Mi 



Y.,Y,k-S%{Mt 



(19) 
(20) 



To compute Zij we thus have to make the connection 
between the 'loop' {hi, /i2, /13) and the 'tree' (/i, H, A) 
states. Neither the zero-search nor the diagonalization 
procedure allow to do this in an unambiguous way. For 
example, level crossings may occur when searching for 
the zeros of det . 

The algorithm currently used by FeynHiggs is: com- 
pute ZHiggs and the associated masses Mi for all per- 
mutations TT of Higgs states involved in the mixing and 
choose the one which minimizes 



M, 



^ I - ZHiggs 



'■rr{i)j\ 



(21) 



where C is the mixing matrix that comes out of the 



Mf^^ , i.e. a by-product 



diagonalization of with q^ 
of the zero-search. 

This is an empirical recipe, so the different dimen- 
sions of AI and Z should not be taken very seriously. 
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r(hi -> T+T-) I MeV 




-77 -Tr/2 7r/2 tt -tt -7r/2 7r/2 



ZHiggs UHiggs(i;2 on-shell) UHiggs(i;2 = 0) 

Fig. 2. The decay width of the lightest Higgs boson to 
r leptons as a function of ^)Xt for tan /9 = 5 (left) and 
tan/3 = 15 (right). The other parameters are: Msusy = 
M3 = Ma = 500 GeV, ^ = 1000 GeV, Afj^± = 150 GeV, 
Xt = 700 e'-^^t GeV. 

The permutation is decided in nearly all cases by the 
mass pattern. The \C — Z\ term becomes relevant only 
for (almost) degenerate masses where it can tell e.g. 
the symmetric from the antisymmetric state. 

The numerical effects of using the various mixing 
matrices in a physical amplitude are shown in Fig. [2l 
For the chosen parameters, UHiggs((7^ on-shell) gives 
results closer to the full result than UHiggs(g^ — 0), 
with deviations at the few-percent level. For a detailed 
discussion see Ref. [B]. 

The mixing of the neutral Higgs bosons in the 
cMSSM has also been added to FeynArts A spe- 
cial version of the MSSM Model File [23, HMix.mod, 
provides two sets of appropriately mixed Higgs bosons: 

- S [0 , {/i>] = UHiggs [/i , i] S [i] , and 

- S[10,{/i>] ZHiggs S[i]. 

The latter is inserted only on external lines. 



3 Benchmark Scenarios 

FeynHiggs has long included Benchmark Scenarios [21 
which are useful in the search for the MSSM Higgs 
bosons. The idea is to vary only Ma and tan /3 and 
keep all other SUSY parameters fixed. 

Constraints such as Cold Dark Matter (CDM) have 
been ignored in these scenarios. It might be desir- 
able to investigate M^-tan /3 planes in agreement with 
CDM and other external constraints, however (if the 
planes are derived in a GUT-based model, see Ref. [22] 
for a discussion). 

For the Constrained MSSM (or mSUGRA) as can- 
didate model, the CDM constraints turn out to be 
too severe, i.e. cut out almost all available parameter 
space. This is different in the NUHM (Non-universal 
Higgs-mass model) [23], where the assumption is that 
there is no unification of scalar fermion and scalar 
Higgs parameters at the GUT scale. As additional free 
parameters in this model one can choose Ma and /i. 

In Ref. [22I four M^i-tan /3 benchmark planes have 
been defined that are in agreement with the CDM and 
other low-energy constraints (see also Ref. fM' ) . From 



a technical point of view, the NUHM introduces non- 
trivial relations between parameters, which thus can- 
not be scanned naively by independent loops. Feyn- 
Higgs 2.6 offers the new format of Parameter Tables 
to deal with such cases. 

Input parameters can either be given in an input 
file (as in previous versions) or interpolated from a 
table, in almost any mixture. The table format is fairly 
straightforward: 



MT 




MSusy 


MAO 


TB 


At 


HUE 


170 


9 


500 


200 


5 


1000 


761 


170 


9 


500 


210 


5 


1000 


753 


170 


9 


500 


200 


6 


1000 


742 


170 


9 


500 


210 


6 


1000 


735 



For two given inputs (typically Ma and tan/3) the 
four neighbouring grid points are searched in the table 
and the other parameters are interpolated from those 
points. An error is returned if the inputs fall outside 
of the table boundaries (i.e. no extrapolation). 

The four predefined NUHM M^-tan/? planes [22] 
can be obtained from www.feynhiggs.de/planes. The 
definition of new planes by the user is possible. 

The Table concept is actually embedded into the 
new FeynHiggs Record. This is a data type which cap- 
tures the entire content of a FeynHiggs parameter file. 
Using a Record, the programmer can process Feyn- 
Higgs parameter files independently of the frontend. 

4 Output of FeynHiggs 2.6 

We give a short overview of the output routines of the 
FeynHiggs library. 

FHHiggsCorr - All Higgs-boson masses and mixings: 
^hi,2,3' ^^H±i "off, UHiggs, ZHiggs. 

FHUncertainties - Uncertainties of the masses and 
mixings. 

FHCouplings - Couplings and Branching Ratios for 
the following Higgs decay channels: 

/ii,2,3 - !],ii,ZZ^*\WW^*\gg, H± ^ /«/', 
h,Z*,hihj,H+H- , hiW^*, 

fifjj fifji 

Y^Y^ Y°Y° Y°Y^ 
A^ A,j : ; AiAj f 

plus the corrsponding channels of an SM Higgs with 
mass Mh,: hf% ^ ff,^^,ZZ^*\WW'-*\gg. 

FHHiggsProd- Higgs production-channel cross-sections 
(SM total cross-sections multiplied with MSSM effec- 
tive couphngs, see Ref. [25] ) 

— gg ^ hi ~ gluon fusion. 

— WW hi, ZZ ^ hi - gauge-boson fusion. 

— W ^ Wh^, Z ^ Zh,- Higgs-strahlung. 

— bb bbhi - bottom Yukawa process. 

— bb bbhi, ~ bottom Yukawa process, one b tagged. 

— ti tihi - top Yukawa process. 
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— tt ^ tthi - stop Yukawa process. 

FHConstraints - Electroweak precision observables, 
see e.g. Ref. [26 for details: 

— Ap at 0{a,aas), including NMFV effects. 

— Mw, sin^ ^eff via SM formula + Ap. 

— BR(6 sj) including NMFV effects [57]. 

— ((7p — 2)susY including full one- and leading/sub- 
leading two- loop SUSY corrections. 

— EDMs of electron (Th), neutron, Hg. 

5 Download and Build 

— Get the FeynHiggs tar file from www.feynliiggs.de. 

— Unpack and configure: 

tar xfz FeyiiHiggs-2 . 6 . 1 . tar . gz 
cd FeynHiggs-2 . 6 . 1 
. / configure 

— "make" builds the Fortran/ C-t- -I- part only, 
"make all" builds also the Mathematica part. 
The build takes about 3 min on a Pentium IV. 

— "make install" installs the package. 

— "make clean" removes unnecessary files. 

The build was tested on Linux, Tru64 Unix, Mac OS, 
Windows (Cygwin) and also with Mathematica 6 (non- 
trivial due to its many incompatibilities) and older ver- 
sions. 



6 Usage 

FeynHiggs has four modes of operation: 

— Library Mode: Invoke the FeynHiggs routines from 
a Fortran or C/C-I--I- program linked with libFH . a . 

— Command-line Mode: Process parameter files in 
FeynHiggs or SLHA format from the shell prompt 
or in scripts with the FeynHiggs stand-alone exe- 
cutable. 

— Web Mode: Interactively choose the parameters at 
the FeynHiggs User Control Center (FHUCC) and 
obtain the results on-line. 

— Mathematica Mode: Access the FeynHiggs routines 
in Mathematica via MathLink with MFeynHiggs. 

All programs and subroutines are documented in man 
pages. 

6.1 Library Mode 

The FeynHiggs library libFH.a is a static Fortran 77 
library. Its global symbols are prefixed with a unique 
identifier to minimize symbol collisions. The library 
contains only subroutines (no functions), so that no 
include files are needed (except for the couplings) and 
the invocation from C/C-I--I- is hassle-free. Detailed 
debugging output can be turned on at run time. All 
routines are described in detail in the API guide and 
on man-pages. 



6.2 Command-line Mode 

The user submits a parameter (text) file, such as 

MT 170 . 9 

MAO 200 

TB 50 

MSusy 975 

Abs(M_2) 332 

Abs (HUE) 980 

Abs(At) -300 

Abs(Ab) 1500 

Abs(M_3) 975 

to the FeynHiggs executable with a command like 

FeynHiggs file Iflags^l 

where the flags are optional. The output is a human- 
readable version of the results. Details of this (rather 
voluminous) output are tagged with a 7, and can thus 
be masked off with 

FeynHiggs file Iflags'i I grep -v 7, 

The table utility converts the output to machine- 
readable format, for example 

FeynHiggs file Iflags'i I table TB MhO > outfile 



The new 'table' statement in the parameter file loads 
the table (see Sect.[3|) and associates two interpolation 
variables with it. The changes are rather minimal: 



Input File "table" 


"inline table" 


MAO 200 


MAO 200 


TB 50 


TB 50 


table file.dat MAO TB 


table - MAO TB 




MAO TB At MUE . . . 




200 5 1000 761 




210 5 1000 753 



Loops over parameter values (parameter scans) are 
possible as in former versions: 



- MAO 200 400 50, linear: 200, 250, 300, 350, 400, 

- TB 5 40 *2, logarithmic: 5, 10, 20, 40, 

- TB 5 50 /6, number of steps: 5, 14, 23, 32, 41, 50. 

6.3 SUSY Les Houches Accord Format 

The FeynHiggs executable can also process files in 
SUSY Les Houches Accord 2 (SLHA2) format [H]- It 
uses the SLHA Library [29] . Processing of SLHA2 files 
can also be done in Library Mode with the subroutine 
FHSetSLHA. 

FeynHiggs in fact tries to read each file in SLHA 
format first and if that fails, falls back to its native 
format. 

6.4 Web Mode 

The FeynHiggs User Control Center (FHUCC) is on- 
line at www.feynhiggs.de/fhucc. It is a Web interface 
for the command-line frontend. The user gets the re- 
sults together with the input file for the command-line 
frontend. A screen-shot is shown in Fig. [3] 
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^0 O B Mozilla Firefox 
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The FeynHiggs User Control Center 


You can still access the version 2.5.I. 
You can siill access the version 232. 




Flags 

Scopeof the 1-loop part: |full WSSM 
I -loop field rcnomialization: | DRbar 
I -loop tan(bcta) rcnorraalization: | DRbar 

Mixing In the neutral Hi££S sector: | 2x2 ChO-HHJ mixing = real parameters _tJ 
Af^jiDximaiion for ihe 14oop result | no approximation _J 




Fig. 3. Screen-shot from www.feynhiggs.de/fhucc. 





6.5 Mathematica Mode 

A more powerful interactive environment is provided 
by the Mathematica interface of FeynHiggs. The Math- 
Link executable MFeynHiggs must first be loaded with 

Install ["MFeynHiggs"] 

and makes all FeynHiggs routines available as Math- 
ematica functions. In combination with the arsenal of 
standard Mathematica functions such as ContourPlot 
and Manipulate, even sophisticated analyses can be 
carried out easily. 



7 Summary: Main New Features 

Version 2.6 of FeynHiggs introduces the following new 
features: 

— Higgs masses are computed as the real part of the 
complex pole. 

— Two kinds of 'mixing' matrices with different prop- 
erties (UHiggs, ZHiggs) are returned. The user can 
choose which mixing matrix to use in all Higgs pro- 
duction and decay channels (default: ZHiggs). 

— Inclusion of the full cMSSM two- loop 0{atas) cor- 
rections in highly optimized form. 

— Inclusion of full one- loop NMFV effects. 

— Possibility to interpolate parameters from data ta- 
bles. Availability of Af^-tan (3 planes in agreement 
with CDM constraints. 

— Total Higgs production cross-sections in effective 
coupling approximation 

— EDMs of electron (Th), neutron, Hg. 
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